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Propagation of electromagnetic waves in multimode waveguides of variable height is investigated. 
The model consists of two uniform rectangular waveguides connected by a linearly tapered waveguide 
section. Using a generalized reciprocity theorem for waveguide junctions, a previous quasi-optic 
solution of this problem is 'extended to account for reflected waves. The results have application to 
the theory of VLF radio propagation when the effective height of reflection of the ionosphere boundary 
varies along the path. The analytical investigation has been complemented by laboratory measure- 
ments taken from a two-dimensional microwave model, and good agreement with calculated results 
was achieved. 



1. Introduction 

The diurnal variation of the electrical properties of the ionosphere results in a change of the 
effective height of the ionosphere lower boundary along the path of propagation. These spatial 
variations of the effect height of reflection modify the field patterns across the waveguide. If 
the change is sufficiently gradual the waveguide theory developed for a uniform earth-ionosphere 
waveguide may be generalized in a fairly straightforward manner [e.g., see Wait, 1964]. Further- 
more, the observed data [Crombie, Allen, and Newman, 1958; Wait, 1959, 1961] are consistent with 
this picture. However, it is now apparent that near sunset and sunrise lines, mode conversion 
effects may be important. It is the purpose of this paper to consider this problem. 

A two-dimensional dual model waveguide described previously [Bahar and Wait, 1964] is the 
basic tool in the present investigation. In this scaled model, rather ideal conditions are assumed. 
At the ionosphere boundary (considered to be sharply bounded), the tangential magnetic field of 
the VLF radio wave is assumed to vanish, corresponding to a reflection coefficient Ri = — 1. The 
earth's boundary is assumed to be perfectly conducting, corresponding to a reflection coefficient 
R g =l. Consideration is restricted to the lower order modes (which account for most of the energy 
of the VLF radio waves at large distances from the source). It has also been assumed in this 
work that 

hi a < < \C 2 \ and (W2) 1/3 ReC > 2, 

where a is the radius of the earth, h is the effective height of the sharply bounded ionosphere, 
k is equal to 27r/free-space wavelength, and C is the cosine of the angle of incidence at the earth's 
boundary. With these restrictions a flat earth approximation can be introduced. Furthermore, 
the effects of the earth's magnetic field have been neglected. 

It has been shown [Wait, 1964; Bahar and Wait, 1964] that, with the assumptions stated above, 
the modal equation governing the propagation of TM modes in the actual earth-ionosphere wave- 
guide is the same as for the TE modes in the dual model waveguide of half -height corresponding 
to the height of the ionosphere in wavelengths. The cosine of the angle of incidence of the pth 
mode on the boundaries of the waveguide, derived from the appropriate modal equation [Wait, 
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1960], is given by 



c P = 



ir(2p-l) 
2kh 



P=l,2,3 



(1) 



where k = 2ir/k is the wave number mentioned above. 

The original curved earth-ionosphere waveguide is represented by the upper (or lower) half 
of the rectangular dual model waveguide (see fig. 1). The upper (or lower) conducting boundary 
of the model waveguide represents the ionosphere boundary and the plane of symmetry (x, z) repre- 
sents the earth boundary. Also TM modes propagating in the earth-ionosphere waveguide are 
represented in the dual model by TE modes. Of course, it is necessary that only waves with 
electric fields that are symmetric about the jc-axis of the model waveguide are initially launched 
into the waveguide; thus the tangential magnetic field at the plane of symmetry (x, z) vanishes, 
which conforms with the vanishing tangential electric field at the earth's surface. It should also 
be noted then that the pth mode in the earth-ionosphere waveguide corresponds to the (2p — l)th 
mode in the model waveguide since the even modes with asymmetric electric fields must be ex- 
cluded. A more detailed discussion of the modeling technique is given elsewhere [Maley and 
Bahar, 1963, 1964; Bahar, 1964]. 

The diurnal change in the effective height of the ionosphere derived from phase velocity 
measurements as reported by Pierce [1955] and Crombie, Allen, and Newman [1958] have been 
discussed by Wait [1959, 1961]. From the viewpoint of mode theory of VLF propagation, Wait 
obtains a value of AA (the change in the effective height) between 16 and 18 km, in good agreement 
with values obtained by Bain et al. [1952], who analyzed the interference pattern of the groundwave 
and the first hop skywave at 16 kc/s. The model waveguide with a half-height ha— 12.7 cm (oper- 
ated at 9 Gc/s) would thus represent the earth-ionosphere waveguide 88 km high (at about 13 kc/s) 
under more or less normal conditions prevailing at night. Similarly the model waveguide with 
half-height h a = 10.16 cm would represent the earth-ionosphere waveguide 70 km high under the 
conditions prevailing during daytime. This paper deals with a special case in which the transition 
between the day and night propagation paths is assumed to take place over a length of several 
wavelengths rather than in the abrupt manner discussed in an earlier paper by the same authors 
[Bahar and Wait, 1964]. The transition region is assumed to be wedge-shaped, and it is bounded 
by uniform waveguide sections of unequal cross sections representing the day and night propagation 
paths. 



2h(x) 
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PERFECTLY CONDUCTING EARTHJ 

ELECTRIC CONDUCTING WALLS 



Figure 1. Cross section of the two-dimensional dual 
model of the idealized earth-ionosphere waveguide 
used in the experimental investigation. 
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2. Formulation of the Problem 



The half-height h{x) of the waveguide, as a function of the distance x along the axis, is ex- 
pressed analytically by 



h{x) = h a -f- [u(x) — u(x — L)] — - — - x 4- u(x — L) (ha — h a ) 



(2) 



where u(x) is the unit step function defined as, 



u{x) z 



fO for x < 



1 for x^O, 



L is the length of the transition path, and h„ and ha are the heights of the uniform waveguide sec- 
tions at the terminals of the transition region. The gradient of the height profile in the transition 
region is restricted in the following analysis by the condition 



* -tan = ^4= < 0.25. 



(3) 



The total variation in height is assumed to be about a wavelength and, thus, the length of the tran- 
sition region is restricted by 



L>4(Arf-A„)«4X. 



(4) 



This is not a serious restriction judging from the nature of the problem under consideration. 

The linear transition described here has been investigated by Wait [1962], who considered a 
TM wave incident on a wedge-shaped transition between two parallel-plate waveguides. These 
solutions are extended in this paper. Here the transition region between the uniform waveguides 
is analyzed as a waveguide junction. 

To derive the scattering of the waveguide modes for the case of propagation in the day-to- 
night path, a quasi-optical method is used in which reflections are totally neglected. This approxi- 
mation is very suitable whenever a low-order mode propagates through a multimode waveguide in 
the direction of increasing cross section, as was illustrated rather rigorously in the investigation 
of mode conversion for the case of propagation of grazing modes in a multimode waveguide with an 
abrupt height discontinuity [Bahar and Wait, 1964]. The above solution is then used to derive 
the scattering when the direction of propagation is reversed (night-to-day path) by invoking the 
reciprocity theorem. In this manner, reflections for the night-to-day propagation path need not 
be neglected. 

In order to apply the reciprocity condition at each discontinuity in the gradient of the height 
profile (x = 0, x — L in fig. 2), it is necessary to derive a generalized interpretation of the reciprocity 
theorem for waveguide junctions. 



FIGURE 2. Wedge-shaped transition region between the 
day and night propagation paths. 
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At the terminal surfaces in the wedge region (ports A, B, C, and D in fig. 2), the tangential fields 
are described by a finite number of unattenuated modes and an infinite number of evanescent 
modes. The electric "basis" field (transverse field pattern) for the nth mode in ports A and B, 
for instance, are denoted by <& A and 4>J respectively. In the problem under consideration 



® A = cos k A y a , | y\ ^ h (l and <$> B = cos v n (p, \ <p | ^ i// 



(5) 



yhere 



» 2h„ 



(6a) 



and 



V n = 



nu 
2t//' 



(6b) 



It will be found very convenient to use matrix notation throughout this work and thus, before 
proceeding with the analysis, the definitions of the symbols for the matrices and their elements 
are stated. 

Let <$> A and <t> B denote the "basis" field row vectors whose elements are $^ and 4>^ respectively, 
and let <I> denote the total "basis" fields row vector, representing all the "basis" fields in both the 
ports A and B, of the junction under consideration. Hence O is defined as 



q>= [<JM<I>*]. 



(7) 



The quantities a A and b A are defined as the nth mode complex wave amplitude (referred to 
port A) of the wave traveling towards and away from the waveguide junction respectively. Thus 
\a A \ is the magnitude of the forward traveling nth mode and arg (a A ) determines the phase of this 
wave relative to the other waveguide modes. The symbols a A and b A are wave-amplitude column 
vectors whose elements are a A and b A respectively. Similarly, with respect to terminal B, a B and 
and b B are wave-amplitude column vectors whose elements are a B and b B respectively. The total 
wave-amplitude column vectors are defined as 



and b 



b A 



b H 



(8) 



Then the electric field at port A expressed in terms of the "basis" fields and wave amplitudes is 

Ej(y, t) = Re {<$> A [a A + b A ] exp (iwt)}. 

Matrices Y A and Y B are diagonal characteristic admittance matrices whose elements Y A and Y B 
are the nth mode characteristic admittances for ports A and B respectively. The total charac- 
teristic wave admittance is defined as 



Y A 
Y B 
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(9) 



The symbols Z A , Z B ', and Z are characteristic impedance matrices. They are equal to the inverse 
of the matrices Y A , Y B , and Y respectively. 

Matrices W A and W H arc diagonal matrices whose elements W A and W B are the power normal- 
ization factors defined by 



W? = ( *" [Q4?dy=ha and W B n = f* [®Wpdtp = p h ty. 
The total power normalization matrix is defined as 



(10) 



w- 



W A 







h a 








W B 







pi4 



(11) 



where / is the unit diagonal matrix. 

The symbol S A is a diagonal matrix whose elements are S A ; S A can be identified as the sine of 
the angle of incidence of the nth mode on the narrow walls of the model waveguide in port A. It 
is given by 



S;! = [1-(A;W]" 2 , 



n= 1, 3, 5 . . . 



(12) 



where k\ is given by (6a). These elements are related to the mode characteristic admittance 
through the equation 



Y A « = YSi, 



(13) 



where Y is the free-space wave admittance. They are also related to the mode propagation con- 
stants {}* through the equation 



« = A«. 



(14) 



where k is the free-space wave number. 

Symbol S AA is a square reflection scattering matrix, whose element S AA is the complex ampli- 
tude of the nth reflected mode when the rath mode of unit amplitude is incident on the junction 
from port A. Similarly, S BB is the reflection scattering matrix related to the waves in port B. 

The symbol S BA is a square transmission scattering matrix, whose element S BA is the complex 
amplitude of the nth mode transmitted through the junction to port B when the mth mode of unit 
amplitude is incident on the junction from port A. Similarly S AB is the transmission scattering 
matrix related to transmission through the junction from port B to port A. 

The total scattering matrix is then defined as 



§AA gAB 

gBA $BB 



(15) 



In terms of the matrices defined above, the scattering matrix equation of the waveguide junction 
becomes 



b = Sa. 

Similarly, the above matrix quantities can be defined corresponding to ports C and D (fig. 2). 
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(16) 



List of Symbols 

a, a A — forward traveling wave amplitude vectors 
a = radius of the earth 

6, b A = backward traveling wave amplitude vector 

h = height of the ionosphere, or half-height of model waveguide 
h n = half-height of "daytime" portion of waveguide 
hd = half-height of "nighttime" portion of waveguide 

k = free-space wave number 
/n, n, p, q, s = integer indices 

u = unit step function 
x, y, z = Cartesian coordinates 

E = electric field 

E= electric field amplitude matrix 

H= magnetic field 

H= magnetic field amplitude matrix 

R = reflection coefficient 

S = scattering matrix 
S AA = reflection scattering matrix 
S BA = transmission scattering matrix 
S A = sine of the angle of incidence for the mth mode 

W— power-normalization matrix 

Y= free-space wave admittance, characteristic admittance matrix 

F= equivalent admittance matrix 

Z = characteristic impedance matrix 

j8 = propagation constant 

€ = dielectric constant 

X = free-space wavelength 

jjl= permeability 
3>, $> A = basis-field row vector 
p, <p, z = cylindrical coordinates. 



3. Day-to-Night Transition Path 

The incident wave traveling in the direction of the axis (from the narrow to the wide uniform 
waveguide as indicated in fig. 2) is assumed to be the TE m , o mode for which the electric field in 
the narrow waveguide is given by 

E z = af n exp (- ikxSfn^Uy) = Q>m exp {- ikxSf n } cos k A n y (17) 

where k A and S A are defined by (6) and (12) respectively. 

For the case considered, the width of the model waveguide increases monotonically (or re- 
mains constant) as the wave advances along the waveguide, and the cross section of the multimode 
waveguide never changes abruptly. Thus, the assumption that reflections may be neglected can 
be applied; this renders S AA = 0, S cc = (actually this assumption was also verified experimentally). 
The waveguide region shown in figure 2 is regarded to possess junctions for the sake of analysis. 
For example: 

Junction I is bounded by the .% = plane and the p = pb cylindrical surface, 
Junction II is bounded by the cylindrical surfaces p = Pb and p = p r , and 
Junction III is bounded by the cylindrical surface p — p c and the plane x = L. 
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Junction II, being a uniform linear wedge region, will transmit the incident waves without 
mode conversion. 

(a) Scattering through Junction I. 

Let x, y, z be the coordinate system connected with port A {x = plane) and let p, <p, z be the 
coordinate system connected with port B (p = pb cylindrical surface). The axes z and z' are par- 
allel as indicated in figure 2. 

The distance Ax (along the x-axis) between the terminal surfaces x= and p — Pb of junction I, 
as a function of the azimuthal angle <£>, is 

Ax = p b (cos <p - cos t//) **"— (i// 2 - <p 2 ) (18) 

where i// is given by (3) and cos i// is replaced by the first two terms of its Maclaurin series expansion. 
The radius pb is related to the parameters of the height profile (2) by the equation 

h (l h Q 

Pb = ~ — T^-p (19) 

sin ifj \p 

with a similar relationship for p c . 

On neglecting reflections, the electric field at p = pt h by substituting (18) into (17), is given by 

E z (p = p b ) = ai exp {-ikS^Ax} 3M « a & exp {-ikp h Sf n VI2} ■ exp {ikpbS*<p 2 l2} cos kfor. (20) 

The electric field in a uniform radial waveguide is expressed in terms of the orthogonal basis 
fields <$> B by the equation 

mi(kp) 

*. = ?« flj|jW « (2D 

where 4>§ is given by (5), and H (2 ] is the Hankel function of the second kind and order v n . The 
form of (21) follows immediately from an earlier analysis of the wedge region [Wait, 1962]. 

It now immediately follows that the continuity condition of the electric field at p — pb is 
given by 

< exp { - ik Pb Sim exp {ikptSi<p*l2}<K - £ Q*b° - £ <*«<. (22) 

n n 

Now, subject to the restriction of (3) on the gradient of the height profile, 

, A miry mir sin <p rmnp 

To solve for the scattering coefficients Sfj*, premultiply by (£) B and integrate with respect to <p 
over the interval (0, i//); thus 

^^^exp{-^^ 2 /2} 



' I exp {ikpbSf n (p 2 /2} cos v n ip cos v m <pd<p 
Jo 



■^exp{-ikp h SiVmif,? n . (24) 
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As indicated before [Wait, 1962], the integral I%£ in the above equation can be expressed in 
terms of a Fresnel integral F(x) defined by 



F(: 



x)=\ X 
Jo 



ITTXT x 

exp i ^ \ ax. 



(25) 



To be explicit, 



7 »"' = i(i 



— i(v m — V n f 



4tt' 



/?<y /2 



*+ 



Vm — Vn 

2ai 



<^)v^ 



+ exp 



— i{v m + V n ) 






??r (*+*&* 



2cd\ 1/2 / v m + v„ 



(26) 



^here 



, kpbSj 



(27) 



(b) Transmission through the uniform radial waveguide Junction II. 

In this junction no scattering takes place, and it can readily be shown [Wait, 1962] that the 
transmission coefficients in the forward and backward directions between the terminal surfaces 
pb and p c are given, respectively, by 



S™ = Hil(kp c )/Hil(kp b )8 n 



(28) 



S™ = Hil(kp b )IHM(k Pc )8 n 



(29) 



The relationship between the coefficients S^f n and S^ checks with the general reciprocity 
theorem (appendix A) applied to uniform radial waveguides. The H^ and Hj, 2) are the Hankel 

functions of the first and second kind (of order v m ) respectively. 

(c) Scattering through Junction III. 

One should appreciate that the wave incident on Junction III comprises several modes even 
when only one mode is incident on Junction I. However, in the analysis it is necessary to con- 
sider only one of these modes incident on Junction III. On using superposition the result for 
several incident modes is readily obtained. The incident wave approaching Junction III from the 
radial waveguide is thus taken to be 



HfHkp c ) 



<£> c 



(30) 



where <£>m = ( &m ls given by (5). 

The basis fields in the rectangular waveguide are 



4>£= cos k D n y 
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(31) 



where 

riTr 



k° = — (32) 



The electric field in the rectangular waveguide expressed in terms of these basis fields is 

£* = £ 6£ exp [-£*(* -L)Sg] (33) 

where 

S£=[l-(*2/*)*P / *- (34) 

Neglecting reflections, and following the analysis applied to Junction I, it readily follows 
that the continuity condition for the electric field at x = L is given by 

ag exp {ikpcSWm exp {-ikp c Sg<p*/2} = 2<^6* = 2*? S S£°£ (35) 

n n 

in which S^ = S^* 

To solve for S££, premultiply by $^ and integrate with respect to y over the interval [0, ha] . 
This gives 

•^nm"" ex P { ik PcS c m Vl2} ^ I ex P {"~ ikpcSmV 2 ! 2 } cos ""^ cos A:£ydy 

~ exp {ikpcSffill}— exp {— ikp c S c m <p 2 l2} cos ^ m </? cos i/„^d^ 
</> Jo 

= exp {tlpeS^/2} | (/»</,)*, (36) 

where, as in (23), it is assumed that 

id a dy dip 

k D n y~v n <p<m& — ~"7 (37) 

and ffi)* is the complex conjugate of /££. The latter is defined by 

/f^ exp {ia£ } ^> 2 } COS V n <P COS Vmifdif, 

Jo 
which in turn can be expressed in terms of the Fresnel integrals as in (26). Thus 



(38) 



«£ = 



2 



(39) 



Combining the results of the scattering and transmission in the separate sections of the wedge 
region and denoting the scattering coefficients of the composite Junction (with terminal planes at 
x = and x = L) 7\ it can be readily shown that the transmission scattering matrix (from port A to 
port D) is given by 

T DA = S DC S CB S BA (4Q) 
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where S BA S CB and S DC are given by (24), (28), and (36). The reflection scattering matrix for port A 
is assumed to be the zero matrix 

T AA = 0. (41) 

The electric field at x = L is therefore given by the matrix multiplication 

E z (x = L) = <& D T DA a A . (42) 

4. Night-to-Day Tapered Transition Path 

When the cross section of the model waveguide decreases as the incident wave propagates 
in a nonuniform waveguide (propagation in the direction of decreasing p or%), the reflections cannot 
be considered negligible unless the gradient of the effective height profile is very small. As in the 
case when the height profile had an abrupt discontinuity [Bahar and Wait, 1964], it was observed 
experimentally that there were substantial reflections in the case of a wedge transition region (for 
the night-to-day transition path), even when the gradient (3) of the height profile was i// — tan i|/ 
= 0.25. In view of this, the method developed in the previous section to derive the scattering 
matrices is not effective in this case. However, since a very good approximation to the solution 
of scattering in the day-to-night propagation path is already available from the preceding analysis, 
the solution to the problem of scattering in the night-to-day path can be obtained readily from it by 
the use of the reciprocity theorem. The reciprocity theorem, for waveguide junctions with uniform 
waveguide ports of constant cross section and a linear axis of propagation, is given by the matrix 
equation 

WYS=WYS (43) 

where the matrices WY and S are defined in section 2 and the curled symbol above the matrices 
represents the transpose operation. Equation (43) states that the matrix product WYS is sym- 
metric. The corresponding relationship between the transmission scattering matrices S AB and 
S BA is 

W A Y A S AB = ^A Y B W B^ (44) 

In terms of the elements of the matrices, the reciprocity theorem is given by 

WA yA CAB = TJVByBQBA ■ //in 

rw m x nt-'mn rr n x n^nnv \*°) 

The derivation of the above theorem is given by Kerns [1961]. A generalization of this the- 
orem for the case when the cross section of the terminal ports of the waveguide is not constant 
is given in appendix A. The scattering matrix S AB can therefore be obtained directly from the 
previous solution for S BA without involving any further approximations, since the reciprocity 
theorem is exact. Hence 

SAB = [JF AYA]-l S BA Y BjyB (46) 

The fields transmitted through the wedge region to port A can be derived in a straightforward 
manner by means of the scattering matrix T AD in the following matrix equation: 

E z (x = 0) = & A T AD a D , (47) 

where the transmission scattering matrix T AD is related to the matrix T DA , given by (40), through 
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the reciprocity theorem as in (46), 

Tad = [W a Y a ]-'T da Y d W d = - Z A T DA Y D (48) 

where the power normalization matrices W A and W D and the impedance and admittance matrices 
Z A , Y A , Y D are defined in section 2. 

In order to derive the fields within the wedge region or to obtain a first-order approximation 
of the reflection scattering matrix T DD , it is necessary to apply the generalized reciprocity theorem 
(appendix A) to each of the subsections of the wedge region, (since ports B and C are not rectangu- 
lar). In particular it is necessary to determine the scattering matrices S AB , S CD , S BB , and S DD . 
This shall be done in the remainder of this section. 

Using the generalized reciprocity theorem, it can be readily shown that the transmission scat- 
tering coefficients through Junction I (fig. 2) from the terminal surface p — Pb to the terminal surface 
x = is given by 



tub yrs CBA 

car— n n CfiA— nm Mm 

m " WAVTA^nm /irkoi \ { '> 

s&(™)l«8GWI" 

where YJ A and YJ B are the total wave admittances defined by (A. 10) and (A. 11) for the uniform 
rectangular waveguide A and the uniform radial waveguide B respectively. The power normaliza- 
tion constants defined by (10) are 

Wi = ha~W* = pi& (50) 



p b « a /ip. 

Assuming the incident field at the terminal p — pb of Junction I is given by <t> B n a B n , the trans- 
mitted field at the plane terminal x = is 

E z (x = 0) = X <KK = 2 n$Z<- (51) 

n n 

The total field (including reflections) at the surface p = pb is approximately 

Eip = Po) = 2 «« e x P {i*St*x) <■ (52) 

n 

This field can be expressed in terms of the incident and reflected modes at the terminal surface 
p — pb by the equation 

E(p= Pb )=2*$8*n+S'Pat (53) 

P 

where S™ is the reflected pth mode due to an incident rath mode at the terminal surface p = pb. 
Hence, the continuity condition on the electric field yields 

2 *?[8„» + S$K ~ 2 *« S »* e *P {*%**}<£. (54) 

V n 
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To determine the reflection coefficients S^, premultiply by <£> B and integrate with respect 
to (p: 

S™ + 8 pm --Jexp {taffi} J J exp { - iafo 2 } cos v p cos v n <pd<p\ Sj* 

= 2|exp{ia^}/^S^ y (55) 



= S( S pn)* S 



AB 
nm 



where {S BA )* is the complex conjugate of S^ which is given by (24). 

Hence, in matrix notation, the reflection scattering matrix is given by 

S BB ~(S BA )*S AB -I. (56) 

The scattering matrix S CD for the waves that are transmitted through Junction III (in the 
direction of decreasing p), and the reflection scattering matrix S DD at the terminal x = L, are derived 
in the same manner as were S AB and S BB , by the use of the reciprocity theorem and the continuity 
condition of the electric field. It can be verified that 

W D Y 70 (rrkn ) 

»w jyc yrc nm n 2 "»* nm ^ ^ 

and 

S DD ~(S DC )*S CD -I (5 8 ) 

where (57) and (58) are analogous to (49) and (56) respectively. 

The reflection scattering matrix T DD for the composite junction (the entire wedge region) can 
now be obtained either by cascading the scattering matrices of the three elementary junctions 
or by direct inspection: 

T DD = S DD + S DC S CB S BB S BC S CD^ ^ 

where the first term represents the reflections due to the discontinuity in the gradient of the height 
profile at x = L, and the second term represents the reflections due to the discontinuity in the 
gradient at x = 0. 

The field at any point within the wedge region may be derived in a similar manner, since all 
the scattering coefficients have now been derived. 



5. Propagation Through a Wedge Region With a Very Small Gradient in the 

Height Profile 

Experimental data obtained from the model waveguide show that the reflection coefficients 
for either direction of propagation can be entirely neglected when the gradient of the height profile 
in the wedge region is restricted by the condition 

ifj « tan iff = d - ^ 0.03. (60) 
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Hence, for the case (A^ — A a ) = 0.76A, L must be about 25X before the reflection coefficients T™ 
can be neglected together with the reflection coefficients T££. 

Provided that the condition on the gradient (60) is satisfied, it is also seen that 



^^ <0 . 3 , 



(61) 



Hence, a further simplification in the computation of the scattering coefficients can be made. The 
exponent in the integral of (24) can be substituted approximately by the first two terms in the 
Maclaurin series expansion such that, for the propagating modes (S* < 1), 



exp {ia^ip 2 } = exp 



ikptSjp 



1 + 



ikSjpbcp 2 



(62) 



The scattering coefficients can now be evaluated in closed form [Wait, 1962J for either direc- 
tion of propagation, since reflections are now neglected in both directions. With the rath mode 
approaching the junction from the narrow waveguide (port A), the scattering coefficients can be 
shown to be given by 



2 



3 7r 2 ra/ 



S»* = exp{-ik Pb Siri2}' 



2 



f ikpffSf^ijj 2 6nm 



exp 



exp 



(-ikhJSffi 
t 2~ 



77" (nr — ir) 
'nrhaSiiiha — ha) (\ 



— (- 1) 2 ; ra ^ n 
2 



2/A 
ihwSUhd— h a ) 32/? ra 



ra = n 



ttIA 



(ra 2 



- rf 



/ i n n-m , 

2 ( — 1) 2 ; mi^ri. 



In the expression for Sfj 4 the following approximation has been introduced: 
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~~2~ 



exp 



ikhaS^J) (\ 



(63) 



since it has been assumed that 



kpbSffi khaSfJj 



< < 1, 



and 



ijj ~ tan i//- 



h d — h a 



From this expression it is clear that the incident mode essentially undergoes a phase shift on 
propagating through Junction I. In the earlier analysis by Wait [1962], this phase factor was 
neglected. 

From the above results it is also seen that, for n^l, the scattering coefficients decrease 

rapidly as — and that very little energy is converted into the higher order modes. It is also 

seen that the magnitude of the scattering coefficients is inversely proportional to the length of 
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the taper region (the scattered power is inversely proportional to the square of the length of the 
wedge region [Solymar, 1959]). 

Examination of (49) and (56) shows that, provided the condition (61) is satisfied, S^ can be 
neglected, since not only are the reflection coefficients S A * negligible, but even the transmission 
coefficients S%£ (n # m) are small in magnitude. 

Finally, it should be pointed out that the results of the analysis in this section should be re- 
stricted to the cases when the incident wave can be assumed to be essentially comprised of lower 
order modes. 



6. Experimental Results 

Calculations have been carried out for the case of the principal mode propagating across the 
night-to-day path with /i a = 3.05A and Ad = 3.81X (fig. 2) for the following castes: 

(a) L = 20k 

(b) L = 10X. 

The transverse electric field variations (amplitude and phase) obtained from the calculated scat- 
tering coefficients are plotted in figures 3 and 4, together with the experimental data from the 
model waveguide. The scattering coefficients T AD are tabulated in tables 1 and 2. 

For the case L = 20\, the approximate formulas (63) for the evaluation of the scattering co- 
efficients between the rectangular and radial sections have proved to be quite accurate. [Equation 
(48) has been used to obtain the scattering coefficient of the entire junction. This involves the 
application of the reciprocity theorem so that reflections are not neglected for the night-to-day 
path.] 

It should be noted that even in this case where the gradient of the height profile is very small, 
the amplitude of the third mode is approximately 12 percent of the principal mode in the narrow 
waveguide section. The amplitude of the higher-order modes m ^ 5 drops substantially (the ampli- 
tude of the fifth mode is about 2 percent of the principal mode). This indicates that for the case 
of nonperfectly conducting boundaries a surface impedance concept would be applicable since 
most of the incident wave is scattered into near grazing modes. 

From the computed data it is obvious that for the case L = 10X the approximate formulas for 
the scattering coefficients (63) are not appropriate and the more accurate formulas (24) involving 
the Fresnel integrals should be used. Special care should be taken in evaluating the transmission 
coefficients in the tapered region (28), (29) since the values of the order and the argument of the 
Hankel functions cross over, with the order being less than the argument for the near grazing 
modes and the argument being less than the order for the higher order modes. In all the above 
calculations only the first six (even) modes have been accounted for (ti=1, 3, 5, 7, 9, 11), this 
being justified by the fact that most of the incident energy is scattered into the lower order modes. 

It is rather revealing to investigate in some more detail the coefficients that are involved in 
the calculation of the amplitude of the third mode transmitted into the day path with the principal 
mode of unit amplitude incident from the night path. This is obtained by the following summation: 

N 
P 

Since it has been noted that only the near grazing modes are of particular significance, in order to 
simplify the following discussion only the first and third mode will be considered; hence 

TAD ~ QABCBCQCD-L. CAB CBCCCZ) 
7 31 ^33^33^31 ^^31 ->11 ->11 • 
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FIGURE 3. Amplitude and phase variations at the trans- 
verse plane x = Ofor the case when a TE 10 mode traveling 
in the direction of the negative \-axis is launched into 
port D with L = 20\ (fig. 2). 



Figure 4. Amplitude and phase variations at the trans- 
verse plane x = Ofor the case when a TE, mode traveling 
in the direction of the negative x-axis is launched into 
port D with L = 10\ (fig. 2). 



Table 1. Scattering coefficients T£, D for the night-to-day propaga- 
tion path /?„ = 3.05\ h d = 3.18k L = 20\ 



TABLE 2. Scattering coefficients T^, D for the night-to-day propaga- 
tion path h a = 3.05X h d =S.Sl\ L= 10A 



n 


Re(Pf) 


Im(Pf) 


1 
3 
5 

7 
9 
11 


0.1051E+01 

.1887E-01 

.7920E-04 

.8953E-02 

-.2212E-02 

-.8585E-03 


0.3583E-00 
-.1351E-00 
.1940E-01 
-.2083E-02 
-.1139E-02 
-.7392E-03 



n 


RefPf) 


lm(T%) 


1 
3 
5 
7 
9 
11 


0.1090E + 01 
-.1194E-00 
-.3798E-01 

.1022E-O2 
- .2000E - 02 

.1086E-O2 


0. 1 765E - 00 

-.1125E-00 

.6171E-01 

-.1364E-01 

.9152E-02 

.5656E-03 



Where E±XY=\0 ±XY 



Where E±XY= 10±- 
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Figure 5. Amplitude variations along the negative x-axis 
for the case when a TE /0 mode is launched into port D 
with L = 10k (fig. 2). 



) -25 -50 -75 -100 -125 

DISTANCE IN WAVELENGTHS FROM EDGE OF DAY PATH X=0 



The first term is clearly the contribution to T^ due to mode conversion at Junction III (fig. 2) while 
the second term is due to mode conversion at Junction I. 

Now it is interesting to note that, as was predicted by the preceding analysis, the amplitude 
of the scattering coefficients S£f and Sgf for L = 20X is half that for the case L = 10X. Nevertheless 
the amplitude of the conversion coefficient T£? is 12 percent of the amplitude of the transmission 
coefficient Tff for the case L = 20A compared with only 13 percent for the case L = 10A. It should 
be noted that this is not due to a large reconversion coefficient from the second mode back into 
the first mode at Junction I for the case L= 10A since |S$f | ~ 1 in both cases where L = 10A and 
L — 20K. The behavior of the coefficient T^f is best understood by noting the phase relationship 
between the two major contributions to T^ (i.e., conversion at Junction I and II). While these 
contributions are only about 10° out of phase for the case L = 20X, they are about 90° out of phase 
forL=10\. 

The above discussion sheds light on design considerations of a taper section (between two 
multimode waveguides of different cross sections), with minimum mode conversion. While it 
has been pointed out in the discussion of the wedge region with the very small gradient that the 
conversion coefficients decrease as 1/L, this is strictly true only if each Junction (fig. 2) is taken 
separately. On treating the composite junction as one unit it is clear that the total conversion 
factor does not decrease monotonically with increasing L. Thus in order to surpress mode con- 
version it is necessary to choose the length of the linear taper L such that the two major contribu- 
tions to T^ are exactly out of phase and thereby tend to cancel out (destructive interference). 
This aspect of the problem is not pursued any further in this paper, which is principally concerned 
with the analysis of the propagation problem. 

The amplitude of the electric field along the axis of the waveguide is plotted in figure 5 for 
the night-to-day propagation path with L = 10\. The undulation of the electric field for x > L is 
clearly seen to be principally due to the interaction of the principal mode with the third and fifth 
modes. 



7. Concluding Remarks 

It is evident that the quasi-optical technique is applicable to problems of propagation in multi- 
mode waveguides when the direction of propagation is from the narrow to the wide waveguide. 
However, for propagation in the opposite direction it is necessary to modify the quasi-optical 
approach. In this paper^ the required extension has been carried out by applying a generalized 
reciprocity theorem for waveguide junctions with nonuniform cross sections. Substantiation for 
this approach is obtained by using a modeling technique which at the same time has provided in- 
sight into the nature of the problem. 

Finally, it might be mentioned that transmission through a transition region with a gradually 
varying height profile (not necessarily wedge-shaped) between two multimode waveguides with 
finite surface impedance boundaries can also be solved by a generalization of the methods described 
here. This case will be discussed in a subsequent paper by the first author. 
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8. Appendix A. The Generalized Reciprocity Theorem for Waveguide Junctions 

The characteristic admittance for a waveguide of nonuniform cross section is generally not 
the same for a forward and a backward propagating wave. It is also a function of the distance along 
the axis of the waveguide. Only in straight cylindrical waveguides with constant cross sections 
is the characteristic admittance the same for both directions of propagation, and independent of 
the distance along the axis of the waveguide. 

For a radial waveguide, for instance, the characteristic admittance for the principal wave 
traveling in the direction of increasing radius p (Junction II, fig. 2) is [Wait, 1962], 

HI WV(kp) 



y + ±_ y v n 

E+ HW(kp) ' ( A - 1 ) 

where Y is the wave admittance in free space and where E+ and H + are the tangential electric and 
magnetic fields, respectively. Here H {2) \u) is the derivative of the Hankel function of the second 

kind and order v n with respect to the argument u. For perfectly conducting walls the argument of 
the Hankel functions is 

nir 

where i// is the azimuthal boundary of the radial waveguide as indicated in figure 2. 

The characteristic wave admittance for the wave traveling in the direction of decreasing 
radius p is, on the other hand, 

r -J = -iF^ (A2) 

r " E- B»{kp) {A - Z} 

where the Hankel function of the first kind replaces the Hankel function of the second kind in (A.l). 
It should be noted that Y+ is equal to the complex conjugate of Y~ , and that for large kp, 

Y+ ^Y~ ^Y 

Now, for the region under consideration, it may be shown that [Kerns, 1949] 

f (E'xH")-(E' f xH')^dS = (A.3) 

where S is the closed surface bounding the region of integration, n is the unit vector normal to S, 
and E\ H' and E'\ H" are any two electromagnetic fields of the same frequency that can exist within 
the source-free region bounded by S. Equation (A.3) is valid in the present context since the elec- 
tromagnetic characteristics of the medium are isotropic. Furthermore, for the present applica- 
tion, S consists of the inner conducting surface So of the waveguide junction plus the terminal 
surfaces Sa and Sb of the two-port. At the terminal surfaces of the junction, the electric and 
magnetic fields can be expressed in terms of the basis fields 0„ (defined in sec. 2), and thus (A.3) 
may be written as 



/. 



X 2 |(1> "' (a '" + b '<»^» {Y t"'» ~ F " 6 " ) ~ 4U "»< + 6 »'> <I> '< (F « «<<- Y n l O \dS = (A.4) 

S A> S B in ii 
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where a' m , b' m and a" m , b" n are elements of the total wave amplitude vectors (8), related to the electro- 
magnetic fields E\ H' and E'\ H'\ respectively. For waveguide ports with boundaries correspond- 
ing to orthogonal surfaces of separable coordinate systems, the basis fields are orthogonal; hence, 
in this case (A. 4) reduces to 



2 (a'm + b 'mWm(Yt n a" n -Y m b'n)-(an+bnW m (Y+ a' m -Y m b' m ) 



: (A.5) 



where W m is the power normalization coefficient defined in (10). In matrix notation this can be 
written as 

(a! + b')W(Y+a"- Y-b") = (a" + b")W(Y + a' - Y~b') (A.6) 

where a is the transpose of a. 

Since the transpose of a scalar is equal to itself, (A.6) reduces to 

d , WY T Sa f = d'^WWSa' (A.7) 

where Y 7 , the total characteristic admittance, is defined as 

yr •= Y + + Y~ 

and b is replaced by Sa, where S is the total scattering matrix defined by (15). 

Now since a and a" can be chosen arbitrarily, the reciprocity theorem for waveguide junctions 
is, from (A.7), 

WY r S = (WY^S) (A.8) 



(WY r )=WY T . 

For completeness, the relationship between the scattering matrix S and the equivalent ad- 
mittance matrix Y of the junction can be shown to be given by 

- r + -y-s f + -f 

The equivalent admittance matrix Y relates the amplitude matrix of the electric field E = a+b to 
the amplitude matrix of the magnetic field H = Y + a— Y~b through the equation H= YE. For the 
case of waveguide ports for which Y + = Y~, (A.8) and (A.9) reduce to the more familiar forms. 
The total characteristic admittance for a radial waveguide is 

4Y 
yi = Y + -\-Y= ^— (\ lm 

n n n irkp m ] n {kp) HM(kp) (A - 1U) 

For waveguides with constant rectangular cross sections, such as port A in section 6 for ex- 
ample, the total characteristic admittance is 

n A = 2Yi = 2YSi (A.11) 

where Sf t is defined by (12). 
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